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Abstract. We demonstrate how a known Whittaker function integral identity arises from 
the t = 0 and q —t 1 limit of the Macdonald polynomial eigenrelation satisfied by Noumi’s 
g-integral operator. 
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The class-one -Whittaker functions are of great interest in representation theory, inte- 
grable systems and number theory (see Gerasimov, Lebedev and Oblezin [8] and references 
therein). Givental [5, 9] gave the following integral representation, which we will take as our 
definition: 
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= TT TT 
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where A = (Ai,..., Ayr) G , x = (xi,..., xyr), X = {xk,i: 1 <i <k < N), XN,i = Xi, and 

N / k k—1 \ N—1 k 

Fx{X) = i X ( X - X ) “ X X . 

k=l Vi=l 2=1 / k=l 2=1 

Whittaker functions satisfy the following two integral identities. 

Proposition 1. Suppose u > 0 and \,v £ with 3?(Aj -|- Vj) > 0 for all 1 < i, j < N. Then 
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The first of these identities is due to Stade [13] , while the second follows readily by appealing 
to the fact that ipxix) = where x' = —XN-i+i- They also follow from [ 12 , Corollaries 3.6 

and 3.7] once one observes that the multiplicative Whittaker functions (x) in [12] are related 
to those defined above via (x) = V’iA(logx), with logx = (logxi,... ,logXAr). 

The Plancherel theory for Whittaker functions (see, for example, [2, Section 4.1.1]) implies 
that for x,y G 


dX^l;_x{x)'^lJ\{y)mNiX) = S{x 


2 /)> 


where the Sklyanin measure mxf^X) is defined as 


^ J —-L 

and 5{x — y) is the Dirac delta function for x = y. Using the above identity one readily observes 
the equivalence of the results of Proposition 1 and the following two eigenrelations. The integral 
operator on the right-hand side of (1) is referred to as a dual Baxter operator by Gerasimov, 
Lebedev and Oblezin [ 6 ], wherein the below result is demonstrated in a rather different manner 
than that followed in this present work. 

Proposition 2. For u > 0 and Im(r(;i) < 0, 1 < z < A", 
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n r(-i^j - iu;j)V’-g(x) 
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The aim of this note is to explain how the eigenrelation in (2) arises as the q ^ 1 limit 
of a certain eigenrelation involving Noumi’s g-integral operator and the Macdonald symmetric 
polynomials. 

Definition 3. The ( 7 -shift operator acts by mapping function f(zi,...,Zi,...,ZN) t 
f{zi, ..., qzi, ..., zxf). For tt G C of suitably small modulus, Noumi’s q-integral operator 91“ 
acts on the space of analytic functions in zi,..., zat as 


^ Yl 


q ^Zi — q izj 


Zi - 


TT- {Izilzj-.q)^^ Y[(T 




2=1 


Here \v\ = vi + ■ ■ ■ + vj^. 

This operator as well as the below eigenrelation it satisfies is due to M. Noumi. It first 
appeared in [4, Proposition 3.24], and further details regarding its derivation are forthcoming 
in [11]. It was subsequently discussed in [3, Section 4] and a proof of the eigenrelation was 
given in an appendix therein by E. Rains. The Macdonald g-difference operators are like¬ 
wise diagonal in the basis of Macdonald symmetric polynomials with eigenvalues of the form 
er[q^^t^~^,q^‘^t^~‘^,. ■ ■ ,q^^), with Cr the elementary symmetric polynomial. While these 
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difference operators must commute with the Noumi g-integral operator there is presently no 
easy way to express them through each other. 

Recall the Macdonald symmetric polynomial P\{z) indexed by partitions A and symmetric 
in the 2 ;-variables with coefficients which are rational functions of two auxiliary parameters 
G [0,1) (see, for example, [2, 10]). We have the following eigenrelation. 

Proposition 4. Noumi’s q-integral operator is diagonal in the basis of Macdonald symmetric 
polynomials so that for any partition A, 



{^iPx){zi,...,ZN)= n 


( 3 ) 


l<i<N 


This eigenrelation can he understood as a formal power series identity in ( or as a convergent 
series identity, provided jCj is suitably small. 


The remainder of this note will be devoted to showing how when t = 0 and g —)• 1 under 
particular scaling this eigenrelation yields that of (2). We will proceed formally. Various uniform 
estimates would be necessary for this particular route to yield a rigorous derivation of ( 2 ). Since 
this identity has already been proved through other means, we do not provide these necessary 
details. 


Our derivation has two steps. The first takes the termwise limit of the Noumi operator, 
yielding the eigenrelation (5). It is this step which we do not fully justify. The second step uses 
residue calculus to rewrite the resulting summation as the claimed contour integral formula. 
From here on set t = 0 and take the following scalings 


Afc = (iV - 2/c + l)e Mog(e ^) + e ^Xk, Zfc = C, =-ue^. 


—e 


q = e 


Furthermore, define the following scaled version of t = 0 Macdonald symmetric polynomial 


N(N-l) N{N-1) 
2 2 


^^^^Px{z), A(e) =-Itt^s ^-e ^log(e/27r). 




= e 


Note, at t = 0, the Macdonald symmetric polynomial has been identified with the class-one 


q-Whittaker functions, as defined in the work of Gerasimov, Lebedev and Oblezin [7] (the term 
“( 7 -Whittaker function” is used to denote different objects by different authors). 


Since we will not be proving a theorem, let us summarize here the result we will formally 
demonstrate. Under the above scaling, assume that a sequence of symmetric function f^{z) have 
a limit as e —)■ 0 to some f{w). Then, formally we will show that 




where the limiting operator is defined in (4) and is identified with the dual Baxter operator 


through Lemma 5. This is the key result of this paper. 

Step 1. Let us consider how (3) scales as e —0. The right-hand side at t = 0 equals 


1 



Observe that from the convergence of the Cg-exponential to the usual exponential (see [ 2 , Sec¬ 
tion 3.1.1]), 
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—ue 


e 


iCq^’’^-,q)oc 
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In [8, Theorem 3.1], it was explained how t = 0 Macdonald polynomials (i.e., g-Whittaker 
functions) converge to Whittaker functions. Certain tail estimates necessary for that argument 
were further provided in [2, Theorem 4.1.7]. The resulting convergence holds that 

Note that in both [8] and [2] there was a mistake - the prefactor for A(e) in the definition 
of V’^(x) should be (as above) whereas in [8, Theorem 3.1] and [2, Theorem 4.1.7] it was 

written as iKjV+ 2 ) ^ briefly explain where this error came from (in reference to the 

proof of [2, Theorem 4.1.7]). Comparing equation (3.8) with (4.25) we find that in (4.25) the 
term A^(^) should actually be A''(x£_|_;^) (here £-hl = N). This error resulted in neglecting £+1 
extra factors of which jives with the difference between the incorrect and correct powers 

(jV_l)(jV+2) _ JV(jV-l) _ ^ 

2 2 — JV . 

Turning to the left-hand side, let us consider how each term in the summation scales: 

q'^^Zi - q'^^Zj _ - vi) + {wj - Wi) 
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gieuii _ 
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{qzijzj; q)ui 1 - 1 - i _ ^-uie+ie{wi-wj) 

1 _ 1 _ 

e{l - i{wi - Wj)) eivi - i{wi - Wj)) 
r(l + i(wj-u;^)) 

r(i + ui + i{wj - Wi)) ^ 

where l.o.t. denotes lower order terms in £. Also note that 


N 


N 
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..., e 
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)=nsh.*( 
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N{N-1) (JV+2)(iV-l) 


i=l 


i=l 


where the effect of Sa,Wi is to shift Wi by a. Putting all of this together we deduce that 


<n“= (tt)MJJ 

j^e(z>o)^ i<j 


Wj -Wi+ i{yj - Ui) -r-r r(l -I- i{wj - Wi)) 


N 


n 


Wj-Wi ^^T{l + Ui + i{wj-Wi)) 

J \ J / / 


n(si,«)''‘ ( 4 ) 


satisfies the eigenrelation 

‘^~'^'4}w{x) = e~'^^ ipwix). 


(5) 


This completes the first step of our derivation. As we already mentioned, various estimates are 
needed in order to turn this into a rigorous proof. 

Step 2. The purpose of the second step is to show how Tl “ can be rewritten in terms of 
contour integrals through residue calculus. We will obtain this as a consequence of the following 
slightly more general result: 

Lemma 5. Let a be a positive real number and consider a symmetric function f defined on the 
set {w G : Im(r/;j) < —a, i = 1,..., N} which is bounded and analytic in each variable in this 
set. Then 


N 


{^-f){w) = 


/ dCsAr( 0 'U^'^™* n r (0 - 

l{a+i^r ijtl 


( 6 ) 
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for u > 0, where sn is the following variant of the Sklyanin measure: 

1 ^ 1 

— 1 
i¥=j 

To see that this implies the desired result, let f{w) = ifwix) and observe that / satisfies the 
necessary boundedness thanks to [2, Lemma 4.1.19]. Making the change of variables ^ i—>■ —i^ 
and combining ( 6 ) with (5) we easily deduce (2). 


Proof of Lemma 5. We will prove the result by expanding the integral on the right-hand side 
of (6) into residues and matching the result with the summation on the right-hand side of (4). 

For i,j = the integrand on the right-hand side of ( 6 ) has singularities coming 

from the factor r(.^j — iwi) at each point of the form fj = iWi — v for v G Z>o (where Z>o 
stands for the set of non-negative integers). By shifting each contour to the left to —oo we 
will pick up each of these singularities, and the result is that the integral will equal the sum 
of the associated residues (The fact that the contours of integration can be deformed in this 
way follows from our assumption on / and the known asymptotics of the Gamma function, 
which give cie“ 2 |Im( 2 ;)|'^ < |r( 2 ;)| < C 2 e“ 2 |Im(z)|’' for Re(z) in some finite interval 

and some ci,C 2 > 0 and r/ G M, see, e.g., [1, (6.1.45)]; observe also that, thanks to the factor 
^ there is no pole at — 00 ). Each of these residues is evaluated at a point of the form 
( 6 , • • •, ■^Af) = (iw^mi - • ■ •, - ^n) for some (mi,..., mN) G {1,..., N}^ and i/j > 0 for 

i = 1 ,..., A", so the right-hand side of ( 6 ) equals 


1 

W\ 


E E 




X r(ir(;mi “ H 


yE* +V^) ^ f{wm^+\Vi,..., Wmt, + iz^Ar ) 


1 






T{iWmi + I'j) 


If nii = mj for some i ^ j then the associated term in the above sum equals 0. In fact, such 
a term has a factor of the form r{i'j — t'j)“^r(r'j — which equals 0 because vj — Vi G Z. 

Hence the above sum is restricted to the case where m* = a{i) for some a G Sn and equals 

■-reSiv i/e(z>o)^ 

N 

n n r(iwz^(j) - Vj - iwzi)) 

i=l 

^ _ 

n r(iu;^(i) - 12* - - Vj - Vi))' 

i<j 


By the symmetry of /, the sum in i' does not depend on u G S’at, and therefore if we choose a 
to be the identity we deduce that the sum equals 


i 2 e(z>o)zv 


1 


f{w + iv) 


n r(i(u;i 

i<j 


n r(i(uzj - Wi) - Vj)) 
i¥=j _ 

Wj) -1^1 + i^j)T{i{wj - Wi) 


Vj + Vi) 


To complete the argument matching the right-hand side of ( 6 ) to the right-hand side of (4) it 
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suffices to show that 


n r(i('Wj - Wi) - Vj)) 

1 __ 

n* n r(i(i«i - Wj) -Vi + Uj)T{i{wj - Wi) - Vj + Ui) 

i<j 


n 

Kj 


Wj -Wi + i{vj - Vi) 


n r(i + 'i{wj — Wi)) 
r(l + Vi + i{wj - Wi))' 


Observe first of all that in the last product on the right-hand side of the above equation, the 
terms i = j equal r(l + Vi)~^ = 1 /( 1 ^^!). Therefore if we write r* = iwi, we need to show that 


TT r(rj -n- Vj)) 
U r{ri - rj - Vi + Vj) 


n 


rj - Vj - Vj + Vj r(l + ri- rj)r(l -k Vj - Vj) 

rj - n r(l -k Vj +ri- rj)r(l + Vi+ rj - n) 


(7) 


Recalling Euler’s reflection formula 


r(i-z)r(z) = ^ 


TT 


sin(7rz) 

and the fact that zT{z) = r(l + z) we have 


( 8 ) 


r(l + ri- r-j)r(l rj - rj) _ -k 

rj - n sin(7r(rj - r*)) 

and 

rj -ri- Vj + Vi 

r(l Vj +ri- rj)r(l + Vi + rj - ri) 

_ sin(7r(rj - rj - Vj + Vj)) r(l - rj + rj + Vj - Vj) r(l rj - rj - Vj + Vj) 

TT T{l + Vj+ri-rj) T{l + Vi + rj-n) 

Using these identities and the fact that sin(7r(a + A:)) = (—1)^ sin(7ra), the right-hand side of (7) 
becomes 


ni-ir-jR 

i<j i^j 


r(l - rj +ri + Vj - Vj) 
r(l -h Vj +ri- rj) 


and hence (7) is equivalent to 


i<j i^j 


r(l -ri + rj + Vi- Vj) r(ri - rj - Vj + Vj) 
T{l + Vj +ri-rj) r{rj-ri-Vj) 


= 1 . 


But using ( 8 ) again the left-hand side equals 


ni-ir-jR 

i<j i^j 


TT 

sin(7r(ri — rj 


_ sin(7r(rj - n - Vj)) 

Vi + Vj)) TT 


with 


(-ir 


K = '^{Vi - Vj) '^{Vi - 2 Vj) 
i<j i^j 

n n n 

= ^ (n - 2 m -h l)vm - (n - 1 ) ^ z/m = 2 ^ (1 - m)vm, 
m=l m=l m=l 


which finishes our derivation since k is even. 
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